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Tactical Alternative Responses (STAR) combat model. It is 
presented as a stand alone program, written in 5IMSCRIFT, 
which can be easily modified as a subroutine for any high 
resolution combat model requiring tactical nuclear effects 
simulation. When presented with a group of targets which 
are deemed suitable for attack by tactical nuclear weapons, 
the program will select units to fire, select proper yields 
for multiple yield weapon systems and assess casualties 
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I. INTRODUCTION 



This thesis presents a nuclear targeting and effects 
program that is intended for inclusion in the Simulation of 
Tactical Alternative Responses (STAS) combat model. It is 
presented as a stand alone program, written in SIMSCRIFT, 
which can be easily modified as a subroutine for any high 
resolution combat model requiring tactical nuclear effects 
simulation. 

When presented with a group of targets which are deemed 
suitable for attack by tactical nuclear weapons, the program 
will select units to fire, select proper yields for multiple 
yield weapon systems and assess casualties among Armor and 
Infantry within the target area. 

A. CAPABILITIES 

The model in this thesis analyzes the given targets for 
destruction by nuclear weapons. In so doing it; 

1. Can handle any number of targets. 

2. Can handle any number of nuclear capable firing 
batteries each of which may have any number of 
different yield rounds available. 

3. Can keep track of any number of friendly maneuver 
units. 

4. Provides for friendly troop safety. 

5. Represents targets as composed only of personnel and 
tanks. 

6. Assesses damage to each individual target element in a 
stochastic manner rather than as an expected value. 

7. Assesses damage caused by blast, heat and radiation. 



B. LI HI TAT IONS 



Seme of the items listed as limitations are such only 
because net all weapons are included in the data base or not 
all possible target elements are included in the targets. 
Host of the following can be input into the model with a 
moderate amount of work. The model does net; 

1. Use free flight rockets, missiles or Air Force 
delivered munitions. 

2. Use strategic weapons. 

3. Consider collateral damage avoidance. 

4. Model any other effects such as NBC warning and 
reporting, optimum time of exit for units caught in 
fallout areas, crossing fallout areas, decontamination 
cr medical evacuation and treatment of casualties. 

C. THISIS FEE VIEW 

Chapter II will briefly review the effects of nuclear 
weapons. The damage mechanisms will be introduced while 
formulas and methodology for building the assessment phase 
of the model will be developed. 

Chapter III will discuss some of the more relevant 
definitions the reader must become familiar with and then 
will shew how the analyst in the field solves the target 
analysis problem. The model will mimic this manual solution 
in order to arrive at the same solution a field anaylst 
would get. 

Chapter IV will discuss the mechanics of how the model 
works without being code specific. Some simplifying 
assumptions will be made and approximations will be offered 
to reduce the processing time. 
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Chapter V will explain the code of the model. Each 
routine and function will be discussed separately. All 
global variables, permanent and temporary entities 'along 
with their attributes, set membership and ownership are 
explained and must be fully understood to comprehend the 
model . 

Chapter 71 will identify possible areas for further 
research. Areas include a method for optimum assignment of 
batteries/yields to targets, standard field operations taken 
during and after a nuclear burst and inclusion cf ether 
elements inxc the target area. 
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II 



E FFEC TS OF NUCLEAR WEAP ONS 



A. INTRODUCTION 

A nuclear explosion, like any conventional explosion, 
results from the very rapid release of a large amount of 
energy in a small space - an energy density. In simple 
mechanics, the liberation of energy is manifested as 
pressure, heat and, in the case of nuclear weapons, ionizing 
radiation. The damage mechanisms of nuclear weapons are 
refered to as blast, thermal radiation and nuclear 

radiation. Damage caused by blast is further divided into 
damage caused by overpressure and by dynamic pressure. 

The heat of the weapon is almost instantly transmitted 
through thermal radiation to every object within 

line-cf-sight constraints. Objects at distances of several 
miles will begin to absorb enough heat to cause combustable 
materials to burst into flame. The pressure of the weapon 
is transmitted through a blast wave which emanates 

spherically from the point of detonation and travels at 
speeds closely related to the speed of sound. Upon 
encountering a target element the blast wave may crush the 
target because of a high overpressure, or it may destroy it 
by translation - the act of tumbling it about on the ground. 

Finally, the nuclear weapon has associated with it the 
harmful, highly penetrating nuclear radiation. The 
explosion emits alpha and beta particles, neutrons and gamma 
rays. Due tc their charge, alpha and beta particles have no 
military significance in weapon detonation and are hereafter 
igno red. 
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Neutrcns are one of the basic building blocks of rhe 
atomic nucleus and have no electronic charge. It is the 
lack of charge which allows deep propagation' and 
penetration. 

Gamma rays are electromagnetic radiation. They are 
identical in composition no light, radio waves, and X-rays 
except fcr wave length and are not easily stopped or 
absorbed . 

Beth neutron and gamma radiation are refered tc as 
ionizing radiation. The biological implication of this is 
that as the radiation passes through living tissue it 
ionizes seme of it, changing the chemical structure into a 
different, non-f uncticnal substance. Essentially, it kills 
cells or inhibits them from functioning as intended. 

B. ENERGY DISTRIBUTION 

The percentages cf total energy appearing as blast, 
thermal radiation, and nuclear radiation depend mcstly on 
the altitude at which the blast occurs. For bursts within 
the lewer atmosphere the percentages are about 50, 35, and 

15 percent respectively. Thus, as with conventional 
weapons, blast is the dominant factor. 

C. BLAST 

Targets are damaged by a crushing mechanism caused by 
overpressure or by a tumbling action caused by dynamic 
pressure . 

1 • O verpressure 

Overpressure is the crushing force applied to a 
target. When an object begins to interact with the blast 
wave, the leading edge of the object is subjected to an 
increase in pressure while the trailing edge is still at 
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ambient pressure. This produces a quick net force away from 
the blast. As the blast wave envelopes the object the 
entire exterior is subjected to an increased pressure . while 
the interior remains at ambient pressure. Thus, a crushing 
force is applied to the exterior in an attempt to damage it. 
Lungs and eardrums are easily damaged by this mechanism at 
overpressures of 6 psi. 

2 . Dyn ami c Pressure 

As the name implies, this is the pressure which is 
associated with high winds. Dynamic pressure can damage 
targets by pushing, tumbling and tearing them apart. 
Targets damaged primarily by dynamic pressure are called 
drag sensitive. With the exception of heavily armored 

vehicles such as tanks, most military materiel is drag 
sensitive. Personnel in the open are drag sensitive and are 
damaged easily by flying debris. 

3. El a st Wave C haract e ristics 

Upon detonation, a blast wave of compressed air 
moves spherically away from the point, of detonation. The 
wave speed is well in excess cf the speed of sound. As the 
wave propagates, the speed diminishes until it exhausts 
itself. The quantification of three blast wave 

characteristics is essential to simulating blast casualties. 
The characteristics are overpressure and dynamic pressure as 
well as the time of arrival of the blast wave. 



a. Overpressure 

Fig . 2. 1 [Ref. 1 : 

distance from ground zero 
overpressure as a parameter, 
weapon hut can be used for 
scaling laws applied. Squatio 
equation which will apply here 



p. 1 13 ] is a graph relating 
to height of burst with 
The figure is for a 1 KT 
any weapon with the proper 
2.1 is the cube root scaling 
and in the next two sections. 
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HEIGHT OF BURST (FEET) 




Figure 2.1 Scaled Values of 1 KT Overpressure. 



An example will 
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Given: an 

Find: The 

Solution: 
is; 



8 KT weapon is detonated at a height of 200 ft. 
peak overpressure at 1200 ft. from ground zero. 
The corresponding height of burst for a 1 KT burst 



h 



1 



h 

w 1/3 



200 

8 1/3 



100 ft. 
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and the ground distance is 



d _ 1200 

1 " w 1/3 ' V 73 



600 ft. 



From Fig. 2.1, at a distance of 600 ft. with a height of 
hurst of 100 ft. the peak overpressure is 30 psi. 

Ref. 1 has two other overpressure graphs. One 
is fcr overpressures of 10,000 to 100 psi and the ether is 
for overpressures cf 15 to 1 psi. 

These problems may be worked in reverse to solve 
for the optimum height of burst if one knows the required 
overpressure to defeat a target. This concept will be quite 
important in chapters III and IV. An example will 
demonstrate the method. 



Given: An 8 KT weapon will be used against a target which 

is defeated by 10 psi overpressure. 

Find: Optimum height of burst. 

Solution: The optimum solution is the one which produces 

the required effect at the farthest distance from the burst. 
From Fig. 2.1 the distance is 1460 ft. and the height of 
burst is 740 ft. From equation 2.1 the actual height of 
burst and actual distance where the effect will be felt can 
be found. 

d = d,*W 1/3 = 1460*8 1/3 = 2920 ft. = 890 meters 



h = h 1 -W 1/3 = 740 *8 1/3 = 1480 ft. = 451 meters 

In this manner, the optimum height of burst for 
this damage mechanism may be found for any target when the 
required overpressure and expected yield are known. This 
will be essential later in defining what is known as the 
governing effect. 
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t. Dynamic Pressure 

Fig. 2.2 [Hef. 1: p. 117] is a graph relating 

distance from ground zero to height of burst with dynamic 
pressure as a parameter. Like Fig. 2.1, it is for a 1 KT 




0 ICO 2C0 300 400 500 600 700 000 900 1,000 1,100 1,200 ,iOG 1,400 

DISTANCE FROM GROUND ZERO (FEET) 

Figure 2.2 Scaled Values of Dynamic Pressure. 



detonation, but the same cute root scaling equation applies. 
Another example will demonstrate the method. 

Given: a 27 KT weapon is detonated at a height of 600 ft. 

Find: The dynamic pressure at 2400 ft. from ground zero. 
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Solution: The corresponding height of burst for a 1KT burst 

is; 



h 



1 




600 

27 



1/3 



200 ft. 



and the ground distance is 



d 



1 



d 

„l/3 

w 



2 4 0 0 
27 



1/3 



800 ft. 



From Fig. 2.2, at a distance of 800 ft. with a height of 
burst of 200 ft. the dynamic pressure is 8 psi. 

In a manner analogous to that offered in the 
section abcve, an optimum height of burst may be fcund for 
targets damaged by dynamic pressure. 



c. Time of Arrival 



Fig. 2.3 [Ref. 1: p. 121] is a 

distance from ground zero to height of burst 
of arrival of the blast wave as a parameter, 
and 2.2, it is for a 1 KT detonation, but the 
scaling equation applies. Again, an 
demonstrate the method. 



graph relating 
with the time 
Like Figs. 2. 1 
same cube rcct 
example will 



Given: a 64 KT weapon is detonated at a height of 2000 ft. 

Find: The time of arrival of the blast wave at a distance 

of 3 miles from ground zero. 

Solution: From previous techniques; 

h = 500 ft. d = 3960 ft. 

From Fig. 2.3 the time for a 1 KT burst is 2.8 sec. and 

1/3 1/3 

t = t 1 *W x/ ' 3 = 2 . 8 • 64 x/ J = 11.2 sec. 
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t.000 2,000 3,000 4.000 5,000 6,000 

DISTANCE FROM GROUND ZERO (FEET) 



Figure 2.3 Scaled Values of 1 KT Blast Have Arrival Time. 



It is interesting tc note that as the yield increases, the 
time of arrival increases which means the speed of the blast 
wave is decreasing. 



D. THERMAL RADIATION 

Thermal radiation is a radiant heat transfer mechanism 
which propagates with the speed of light. Upon burst every 
object with line-of- sigh t tc the burst begins to absorb heat 
an a rate predicted by the thermodynamic laws governing 
thermal radiation. Ccmbustable materials will ignite if the 
rate and exposure time are sufficiently high. Personnel 
will suffer first, second, and even third degree burns if 
reaction times are slew and absorption rate is high. Fig. 
2.4 [Ref. 1: p. 291] shows the relationship between yield, 

distance from ground zero and the expected heat absorption. 
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Figure 2.4 Heat Absorption as a Function of Range and Yieia. 



0 



1561 



Yield 

Cal 



.5051 

.4385 



( 2 . 2 ) 



Sguation 
Nuclear 
radiation . 
destructive 
light which 
conditions. 



2.2 is the linear regression for Fig. 2.4. 
bursts emit two distinct pulses of thermal 
The second one contains aoout 99 percent of the 
power. It consists of visible and infrared 
is poorly attenuated by normal atmospheric 
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Any condition which affects the visibility or 

transparancy of the air will modify the transmission of 
thermal radiation. Clouds, smoke, fog, rain and snow will 
absorb and scatter the energy. Also, a cloud cover above 
the burst can reflect energy back to the target and increase 
the therial radiation that would have otherwise traveled 
harmessly into the sky. Obviously, any solid object which 
will not become a victim of thermal radiation will provide 
adequate protection. Such objects include, hills, fcxhcles, 
bunkers, vehicles, trees and other people. 

As stated earlier, the exposure time would be a factor 
in determining total exposure. However, in actual tests it 
was discovered that ccmbustable materials such as wood frame 
houses charred almost instantaneously. Also, due to the 
short duration of the thermal pulse, a steady state thermal 
heat transfer condition was not achieved and the buildings 
were able tc absorb and dissapate the initial burst of 
energy and avoid ignition. 

The point of the preceeding paragraph is that reaction 
time of a potential absorber, such as a soldier on the 
battlefield, will play an insignificant part in determining 
thermal radiation effects. 

Fig. 2.5 [Ref. 1 : p. 314] shows the futility of avoiding 
burns if one is sufficiently close to ground zero. The 
following quote [Ref. 1: p. 313] sums if up; 

At the lcwer energy yields the thermal radiation is 
emitted in such a snort time that no evasive action is 
possible. At the higher yields, however, exposure tc 
much cf the thermaj. radiation could be avoided if 
evasive action were taken within a fraction of a seccnd 
of the explosion time. It must be remembered, of 
course, that even during this short period a very 
considerable amount of thermal energy will have been 
emitted from an explosion of high yield. 
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TIME ISECONOS) 

Figure 2.5 Energy Emmission as a Function of Time. 

E. HOCLZiB BADIAT ION 

1 . Neu tro n Rad iati on 

Dpon detonation cf a nuclear weapon a large quantity 
of free neutrons is released. The neutron release process 
is completed in less than a millionth of a second (Hef. 1: 
p. 340]. This is where the term prompt radiation comes 

from. However, due to collisions in the early stages cf the 
detonation, actual escape may be delayed for a thousandth of 
a second - an insignificant comfort to potential targets. 

Since the neutron has no charge it does not have a 
direct ionizing effect. The biologically harmful ionization 
occurs when the neutron collides with hydrogen in body 
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tissue. The neutron transfers its energy to the nucleus and 
frees it from the accompanying electron. Thus a free 

protron and a beta particle are allowed to cause ionizing 
tissue damage. Previously, beta had been disregarded as a 
source of concern. This was due to the impossibility of 
getting the particle into the body. In the context of the 
current discussion, it is already inside the body. 

Fig. 2.6 [Ref. 1: p. 346] is a graph relating yield 
to the slant range tc the burse with radiation absorption in 




2 5 10 20 50 100 

EXPLOSION rtELO <KT) 



Figure 2.6 Neutron Radiation. 

body tissue as a parameter. For a given yield and slant 
range, the absorbed radiation can be predicted. Eguation 
2.3 is the regression equation for Fig. 2.6. 



N 



5 . 896 



2368.47 
(R-15 0 • In (Yield ) 



(2.3) 
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2 . Gamma Ra dia t icn 

Upon detonation of a nuclear weapon, gamma radiation 
is emitted. There are several mechanisms which cause gamma 
radiation at various times during and after the detonation. 
Ref. 1 (Sect 8.08 - 8.19) contains a detailed explanation of 
the mechanisms. Rather than delve into the sources of gamma 
radiation it is more important to predict absorption rates 
and figure out how to avoid absorption. 




exptosiCN ’x") 



Figure 2.7 Gamma Radiation. 



Fig. 2. 7 [ Ref. 1: p. 333] 
For a given yield and slant range, 
gamma radiation can be predicted, 
regression equation for Fig. 2.7. 



is similar to Fig. 2.6. 
the absorbed dose due to 
Equation 2.4 is the 
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G = 46166 



. 5599 • /Yield . 00 01 ‘Yield • R 

e • e 

. 0053 -R 



(2.4) 



Since gamma radiation is just light in a certain 
range cf the electromagnetic spectrum, it behaves in 
accordance with the same laws, such as absorption. However, 
commcn opaque objects such as steel, concrete, paper, etc. 
which reflect or totally absorb visible light, only absorb a 
portion of the gamma radiation. A common example of this is 
traveling into a tunnel with the radio on in the car. The 
reception weakens as the amount of shielding increases. 
Table I lists seme common materials which could be used as 
protection. 



TABLE I 

Common Gamma Absorbers 



Material 


Tenth Value 
Thickness (inches) 


Steel 


3.3 


Concrete 


11 


Eart h 


16 


Water 


24 


Weed 


40 



i 1 



The Tenth Value Thickness is that thickness of the 
specified material which will cut the radiation by a factor 
cf 10. Thus, 3.3 inches of steel or 16 inches of earth will 
reduce a 500 rad. cose to 5 0 rads. The addition cf another 
3.3 inches cf steel will further cut the dose to 5 rads. 
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Thus , 



D = - D 0 *10' T/Tl 

1 0 T/T 1 



where ; 

D is the predicted absorbed dose 
E 0 is the unprotected dose 

T is the thickness cf protective material. 
Tj is the tenth value thickness. 
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III. N0CL5 AB TARGET ANALYSIS 



A. DEFINITIONS 

Eefore beginning to make a weapon selection it is 
necessary to understand the terms peculiar to the subject. 
Some of the more important ones are presented here while the 
more cbvicus ones are found in the glossary. 

1 • Rad ius oj Tar get (RT ) 

This is the actual radius of the target or an 
equivalent radius if the target is equatable to a circle. 
Ref 2 contains a nomograph to determine the appropriate RT 
for targets equatable to a circle. More importantly, the 
target elements within the target are assumed tc be 
uniformly distributed. 

2. Ra d ius of Da m age (RD ) 

This is the distance from ground zero at which a 
single target element has a 50% chance of receiving the 
specified degree of damage. The RD is a function cf weapon 
yield, height of burst, casualty criterion and, in the case 
cf personnel, the protection level. Fig. 3.1 shews a circle 
drawn with a radius from ground zero such that as many 
damaged elements (dark dots) are outside the circle as there 
are undamaged (light dots) inside the circle. Thus, an 
element at RD is said to have a 50% chance of becoming 
damaged. 
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3 ♦ Minimum Ra d i u s of Damage ( MIN RD) 

If the firing of a large number of rounds against a 
target was simulated, a probability distribution for the RD 
would begin to emerge. Fig. 3.2 is a cumulative 
distribution function for a hypothetical RD distribution. 
The MIN RD is as marked and is equalled or exceeded 90% of 
the time. 




Figure 3. 2 CDF of RD. 
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• C ir e ula r Error Probable (CEP) 

If a large number of rounds obeying a circular 
normal distribution were fired at ground zero and the impact 
points plotted, then CEP would be defined as that distance 
from ground zero that is exceeded as often as not. Fig. 3.3 
shows an idealized plot of such firings. A circle of 1 CEP 
includes 50? of the rounds. 




Figure 3.3 Circular Sormal Distribution of Impact. 

The distribution of impacts is assumed to have a 
bi-variate normal distribution about the aim point with 
equal variances and a zero correlation coefficient. Thus, a 
relationship between CEP and variance can be determined and 
is; 



CEP 



1.1774*a 



(3.1) 
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5. C ir cular Distribution 90 (CD9Q) 

CD90 is very much like CEP except that 90% of the 
rounds are expected to fall within 1 CD90 of ground zero. 
When a priority target is being considered it is CD90 which 
is used as the dispersion parameter. The relationship 
between CD90 and CEP is 



CD90 = 1.83 * CEP 



(3.2) 



6. Pro bab le Error H eiq ht-o f -bu rst (PEH) 

Similar to CEP, but in one dimension, the PEH is 
that distance above and below the desired height of burst in 
which 50% cf the rounds are expected to function. Vertical 
delivery is assumed to be normally distributed abcut the 
desired height of burst. Therefore, the relationship 
between PEH and standard deviation is easily calculated as; 



a 



PEH 

.67 



(3.3) 



Most fuzing mechanisms are highly reliable and accurate. 
This ensures optimum height of burst and maximum casualties. 
It also allows seme simplifying assumptions in Chapter IV. 

7 . Gov ern ing E ffect 

Most targets are affected by more than one damage 
mechanism. In order to simplify the analysis the range at 
which a target is damaged is tabulated and the effect 
corresponding to the largest range is selected as the 
governing effect. Table II lists hypothetical ranges for 
effects for exposed personnel and personnel in open 
foxholes. 
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TABLE II 

Governing effect 



Effect 




Range 


(met ers) 






Ex posed 


Fox hoi- 


Overpressure 


575 


575 


Dynamic 


pressure 


225 


N/A 


Thermal 


radiation 


950 


N/A 


Nucle ar 


radiat icn 


600 


250 



For exposed personnel, it would appear that thermal 
radiation is the governing effect. However, due tc the 

unpredict ibility of target posture, thermal radiation is 

never used as the governing effect [ Ref . 3: P. 23]. Also, 
it is almost universally true that the effects of 
overpressure are felt at a greater range than the effects of 
dynamic pressure. This leaves only two criteria to compete 
for governing effect - blast due to overpressure and nuclear 
radiation. In table II, the governing effect for exposed 
personnel is nuclear radiation and for personnel in foxholes 
it is overpressure. 

8. Desired Height of Burst (DH03) 

The DHOB is the maximum of two possible HOB. 

It is usually desirable to obtain a low air burst 
without causing fallout. Eg. 3.4 is the formula for 

calculating the height of burst fallout safe (HOBfs) for 
yields less than 100 KT and Eg. 3.5 is the formula for 
calculating the HOBfs for yields greater than 100 KT. 
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HOBfs = 3 0 * W 



(3.4) 



1/3 



HOBfs = 55*W 1/3 (3.5) 

W is the expected yield plus 10%. 

Obviously this will produce fallout in 50 % of the 
rounds. Therefore, a safety buffer of 3.5 PSH is added to 
HOBfs. This quantity is now refered to as H0399 - the HOB 
at which there is a 99% chance cf no fallout. 

Depending upon the type cf target and command 
guidance, there is a HOB which will cause optimum damage 
(HOBopt) . This was seen earlier in the discussion on 
overpressure and dynamic pressure and the governing effect 
section above. The desired HOB (DH03) is the maximum of 
HOB99 and HCEopt. In most cases, HOBopt is greater than 
HOB99 and is, therefore, the listed HOB. It is only when 
?EH becomes large that HOB99 is greater than HOBopt and the 
DHOB must be raised above HOBopt, thus reducing expected 
coverage . 

9 . Ground Z ero (G2) 

The point on the ground directly below the position 
where the round actually detonates. 

10. Cs sir ed Ground Zero (DGZ) 

The actual aim point. The DGZ might not be the 
center of the target as will be seen later. GZ is only a 
single realization of the DGZ and is not known until after 
the detonation. 
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B. EAMAGE ESTIMATION 



1 . Coverage 

To cause maximum damage, rhe round must be targeted 
at the cenrer of the target. For various reasons, it may be 
necessary tc deliberately shoot at a location other than the 
target center. If RT and RD are of rhe same order of 
magnitude, then a less in coverage will result as shown in 
Fig. 3.4. 





A) Less cf coverage B) No loss of coverage 



Figure 3.4 Loss of Coverage Due to Offset. 
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The second figure is an expected coverage. There is a 50® 

chance of obtaining the listed coverage of the specified 
damage when using this figure. 

2 • Cheesing the Correc t Sound 
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system and the targ 
determined by target 
weapon yield. Fig. 
in the open, latent 
yield . 



et radius. The appropriate table is 
type, desired effects, weapon type and 
3.5 is a typical table for personnel 
lethatities, short range cannon, 1 KT 



COVERAGE TABLE SHORT RANGE CANNON 

(Distances In meters) 1 VT 



S E D PERSONNEL • LATENT LETHALITY 
LOW AIRBURST 

EFFECTIVENESS 
RADIUS OF TARGET 

TOO 800 900 



EXPO 



RANGE 


500 


600 


2000 


. 96/. 98 


. 93/. 96 


3000 


96/. 98 


. 91/96 


4000 


. 36/. 98 


. 91/. 95 


5000 


. 96/. 98 


. 91/. 95 


6000 


. 96/. 98 


. 91/. 95 


7000 


. 95/. 97 


. 91/. 94 


8000 


. 96/. 96 


.86/. 92 


9000 


. 92 /. 96 


. 78/. 88 


10000 


38/94 


. 71/. 85 



. 85 /. 


.88 


. 31/55 


. 41/. 43 


. 85 /. 


,37 


. 52/. 56 


. 41/. 43 


. 85 /. 


,87 


. 53/. 56 


. 41/. 43 


84 / . 


.87 


. 53/. 57 


. 41/. 43 


. 81 / 


86 


. 53/. 57 


. 41/. 43 


. 80 /, 


.85 


. 55/. 57 


. 41/. 43 


. 75 /, 


,83 


. 54/. 57 


. 41/. 43 


. 68 /, 


.78 


. 53/. 56 


. 37/. 41 


. 61 /, 


.75 


. 50/. 55 


. 36/40 



. 29/32 . 21/. 22 . 09 /. J 2 

. 30/. 32 . 20/. 22 . 10/. 12 

. 30/. 32 . 20/. 22 . 10/. 12 

. 30/. 32 . 20/. 22 . 10/. 12 

. 30/. 32 . 20/. 22 . 10/. 12 

. 30/. 32 . 20/. 22 . 10/. 12 

. 30/. 32 . 20/. 22 . 10/. 12 

. 28 / .32 . 20/, 22 . 08 /. II 

. 26/. 31 . 18/. 21 . 06/. 11 



PROB 

MIN 


EXPT 




RD 


RD 


CD 


460 


526 


80 


458 


532 


93 


454 


532 


99 


447 


538 


111 


448 


540 


128 


442 


540 


141 


434 


533 


176 


420 


520 


205 


400 


490 


234 



ACCURACY DATA 



CEP 


HOB 


PEH 


38 


49 


9 


44 


59 


12 


47 


59 


12 


54 


69 


14 


62 


82 


17 


73 


90 


19 


86 


101 


22 


too 


U 6 


26 


115 


132 


30 



Figure 3.5 Exposed Personnel, Latent Lethalities, SRC, 1 KT. 



Entry arguments are the range (rounded up to the 
next highest listed value) and target radius. If the target 
radius is less than the first listed radius, then the first 
listed coverage is taken as the expected coverage. If the 
target radius is larger than the last listed radius then the 
weapon system is considered unsatisfactory. Otherwise, 
interpolation is performed on the radius and the resulting 
coverage is reported as the expected initial coverage. 

EXAMPLE. If a short range cannon weapon was employed 
against personnel in the open with a target radius of 690 
meters and the range to the target was 5200 meters the 
expected coverage would be 86.9%. 
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2 . 



Numerical Method 



When it is necessary to deliberately shoct -at a 
target with an aim point other than the target center, the 
numerical method must be used. Fig. 3.6 was the chart 
designed for usage with expected coverage, but has been 
ruled invalid by USANCA. (Fig. 3.6 may still be used for 
expected coverage with no offset.) Therefore, when 

employing an offset aimpcint it is necessary to use the high 
assurance chart as shown in Fig. 3.7 which is used in the 
following discussion. 

The numerical method has the following 6 steps. 



1. Calculate RD divided by RT and locate on the left axis. 

2. Calculate CD90 divided by RT and locate on the bottom 
axis. 

3. The intersection of the values predicts the coverage 
withcut offset. 

4. Calculate d divided by CD90 and locate on the left 
axis. 

5. Measure the distance from the point located in 4 above 
to the line labeled "displaced DGZ a/CD90". This is 
refered to as the "measured distance". 

6. Apply the " measured distance" to the riqht of the point 
found in 3 abcve. This final point predicts the 
coverage due tc offset. It must also be remembered 
that this coverage has a 90% chance of occurring since 
the scluticn was found by using the high assurance 
table. 



GIVEN: 



The 


followin 


g example 


is dep: 


.cted 


in Fig. 3.7. 


RD = 


600 


met 


ers 


CD90 = 


100 


meters 


ET = 


800 


met 


ers 


d 


300 


mete rs 
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Find: The fractional coverage with the DGZ displaced. 

Solution: 



Step 


1 


BD/RI 


= .75 


Step 


2 


CD90/RT 


= .125 


Step 


3 


Initial 


coverage = .57 


Step 


4 


d/CD90 


= 3.0 


Step 


5 


Me a sure 


distance 


Step 


6 


Final 


coverage - .48 
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Figure 3.7 



CD90 





IV. THE MODEL 



A. IHTBCEUCTION 

Before the medal is invoiced, a serious condition must 
exist on the battlefield - one in which the commander with 
responsibility and authority for employment of nuclear 
weapons has decided tc use them. This decision must be mads 
externally tc the model. Once the decision is made, the 
model will make all of the decisions necessary to solve the 
target analysis problem, deliver the rounds and assess 
damages to target elements. 



E. I SPOT 

As with any military operation, there are several 
essential elements of information which must be available. 

1. The location and size of friendly maneuver elements 
must be known. 

2. The location and type of nuclear capable artillery 
batteries along with their nuclear load must be known. 

3. The location size and type of enemy targets must be 
estimated. 

4. The commander's guidance for minimum acceptable damage 
to enemy units must be announced. 



C. STEPS OF THE MODE! 

Step 1. For each target select feasible battery/yield 
combinations. 

A. Use the index method first. 
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E. 

C. 

D. 

Step 2. 
Step 3. 

A. 



If the index method provides adequate coverage, 
check for troop safety. If offset is cot 
necessary, file as a solution. 

If offset is required, use numerical method. If 
coverage is adequate, file as a solution. 

Of all feasible solutions, select the one which 
provides maximum coverage. 

For each target , employ the best round. 

Assess damage for each target. 

Assess damage to each tank. 

1. Determine the distance from G2 where a tank 

has a 50 % chance of damage due to 

o verpressure . 

2. Declare each tank as alive or dead as a 
result of a Monte Carlo draw. 

Assess damage for each troop. 

1. Determine the distances from GZ where a troop 
has a 50 & chance of survival due to thermal 
radiation, overpressure and dynamic pressure. 

2. For each troop in the target, perform a Monte 
Carlo draw for each of the above damage 
mechanisms. If any one fails, kill the 
troop . 

3. For each troop in the zarget which survives 
tc this point, determine nuclear radiation 
absorption. Place the troop in one of the 6 
categories of; 

a. Dead 

b. Permanantely incapacitated. 

c. Temporarily incapacitated. 

d. Impaired. 

e. Latent impaired. 

f. Safe 
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Tc dc all of this the model uses various routines and 
functions which are detailed in chapter V. The remainder of 
this chapter will explain how the model proceeds, what 
assumptions are made and, in a few cases, where seme 
simplifying approximations are made. 

D. SELECT EEASIELE E ATT EHY/YIELD COHBIHATIOHS 

The model loops through every target with every possible 
fcatt er y/yie Id combination in an attempt to compile a listing 
of adeguate solutions. The proceedure is tc; 

1. Use the index method. 

2. Determine if offset is necessary. 

3. Use numerical method if offset is necessary 

4. If coverage is adequate, file as a feasible solution. 

1 • Index Method 

Tc start, the range to the target is computed. If 
it is farther than the maximum range of the weapon system 
currently being looped through, then it cannot be engaged by 
that battery. The model will new loop to the next battery. 

Technically, a target which is out of range may 
still be attacked successfully. In this situation an offset 
distance equal to the difference between the actual range 
and the maximum weapon range is used in a numerical method 
analysis. The model does not allow for such an occurrence. 

If the target is within range, the model loops 
through each yield of the weapon system. Starting with the 
smallest yield, the target radius is compared to the maximum 
table target radius. If if is too large, then the model 
loops to the next higher yield. If it is not too large, it 
performs an index method analysis. 



44 



Should the index method analysis provide a coverage 
which is egual to or greater than the command guidance for 
the target type, then a potential solution has been found. 
If the coverage is inadequate, then the model loops to the 
next higher yield. 

It must be remembered that the target/battery/yield 
sequence is a nested loop and it may not always be possible 
to lcop to the next higher yield. At a time when the 
highest yield for a weapon is being looped over, looping to 
the next higher yield actually means looping to the next 
battery and starting again with the lowest yield. The same 
applies tc looping tc the next battery when the model is at 
the last battery. 

Whenever a feasible solution is found with the index 
method, the target is immediately checked fcr necessary 
offset. Thus the loop in temporairly interrupted in order 
to check for possible offset. 

2. C-sJuJ-iiLH De sired Ground Ze ro 

Hecall that the index method required the desired GZ 
to be the center of the target. Any offset necessitated by 
troop safety considerations would require the use of the 
numerical method which would analyze the target as a high 
priority target and result in a degraded coverage. In the 
case where initial coverage was adequate, further analysis 
must be dene to see if offset is required and what effect on 
feasability it has. 

Tc perform an offset is easily done on a map with a 
compass, but is not as easily done with a digital computer. 
Fig. 4.1 shews how it is done manually. The method is to 
draw circles about each friendly maneuver unit center. The 
radius of each circle is the unit radius plus the minimum 
separation distance (MSD) to insure troop safety. The 

desired ground zero must lie outside each circle. The 
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Desired GZ 




© © © 

Friendly Units 

Figure 4.1 Manual Solution to Offset Problem. 
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chosen aim point is the point on the map which is closest to 
the target center and still exterior to every circle. - 

Simply stated, one just gets as close to the target 
center as possible while staying far enough frcm each 
friendly element to ensure safety. To do this envclves 
solving a ncn-linear programming problem with a non-ccnvex 
feasability region. The problem formulation is; 

MIN (X-X m ) 2 + (Y-Y ) 2 
T T 

st. (X-X.) 2 + (Y-Y.) 2 > (MSD + Unit Radius.) 2 

11 i 

i = 1, • • • , No. of Units 

Rather than solve that problem, it is easier to 
solve an iterative geometric approximation. In many cases 
the apprcxima tica is exact. The approximation is as 
follows . 

1. Find the friendly unit whose perimeter is closest to 
the target center. If it is within MSD of the target 
center, move the desired ground zero directly away 
frcm it an amount necessary to satisfy safety. 

2. Repeat the process until; 

A. No unit lies within MSD or 
E. 10 steps have been tried. 

If a solution cannot be found in 10 steps it is 
probably because no feasible solution exists. Fig. 4.2 is a 
typical example. Cf the five units, three of them lie 
within MSD ana the desired ground zero must be moved. Unit 
4 is closest and the desired ground zero is moved directly 
away from it. The new desired ground zero is acceptable. 
It is also exact and not an approximation. This is the 
typical case. 
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Figure 4.2 A Geometric Solution. 
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If unit 2 moves forward as shown in Fig. 4.3, the 
solution requires two iterations and is now an 
approximation. The first iteration moves the desired ground 
zero tc position 1. However, unit 2 is within MSE and 
desired ground zero must be moved again. This time pc-sirion 
2 is selected. It is feasible and close to the optimum 
location . 



3 . Num eric al Method 

The numerical method was discussed in Chapter III, 
Section C-2. This method was developed for use by a nuclear 
target analyst in a field environment using tools commonly 
available in a Divisional Fire Support Element. Such tools 
are F M 101-31-3 [Bef. 2], a compass, dividers and a pencil. 
Thus, the transposition of the "measured distance" and 
reading of a point on a graph seem like simple tasks and 
natural things tc do manually. The method followed in the 
computer code is a slight modification of the manual method. 

The scales in Fig. 3.7 are logrithmic. Therefore, 
the addition of the "measured distance" is really a 
multiplication operation such as; 

CD 90 /RT + — CD90/RT * f(d/CD90) 

where f(d/CE90) is the appropriate multiplier and is only a 
function of d/CD90. It is much more reliable in the field 
to allow the target analyst to use dividers than to require 
him tc correctly multiply. 

In the example on page 39, the initial CD90/RT was 
.125 and the displaced CD90/RT was .450. Obviously the 
appropriate multiplier was 3.60. 

Fig. 3.7 has been included here as Fig. 4.4 with 
point number 7 added. If the value of the point number 7 is 
multiplied by 100, it will be 3.60. Thus, the appropriate 
multiplier can be found directly from the figure without any 
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Position 2 





Figure 4.3 Offset as an Approximation. 
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further modification. Several points on the d/CD90 scale 
and the corresponding f(d/CD90) are in the data base. For 
any d/CE90 the correct f(d/CD90) can be found with one 
interpolation. 

The coverages from Fig. 3.7 are input into the data 
base as a 2-dimensional table. Thus with RD/RT and CD90/RT 
known, a double interpolation will produce the predicted 
coverage. 

4 . Choosing the Eest 3a ttery/Y ield Combina tion 

If, after all of this work, the coverage is now 
inadequate, the solution is abandoned and the program moves 
on. However, if it is adequate, a solution containing the 
target, battery, yield, desired 3Z and predicted coverage 
can be generated and filed away for later use. It would be 
useful to file the feasible solution in such a way that 

recovery of the top one on the list is the best solution. 

Table III shews hypothetical results of an analysis 
requiring offset for the rounds of the MRC battery. Due to 
safety imposed offset, the largest yield, while still 

meeting command guidance, is clearly not the best choice. 

In fact, neither of the rounds from the MRC battery is the 
best choice. The SRC battery with a 1 KT yield will provide 
the highest coverage. Therefore, this program will file 
solutions in descending order of coverage. Should a tie 
occur, the SRC takes precedence over the MRC for 
conservation of force reasons. 



E. EMPLOYING THE ROUND 



When all of the targets have been analyzed for each 
battery, the target analysis portion of the model is 
finished. It is time to select the best battery/yi eld 
combination, employ the weapon and assess damage. Actual 
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TABLE III 

Degraded Coverage Due to Offset 



Eatter y 
Type 

Short Range 
Cannon 



Kedium 

Cannon 



Hangs 







Initial 


Offset 






coverage 


coverage 


.2 


Kt 


41% 


4 1% 


1 


KT 


78% 


78 % 


2 


KT 


82% 


73% 


8 


KT 


97% 


44% 



burst parameters are normally distributed. The actual 
ground zero is picked from a bivariate normal distribution. 

GZ = N (desired ground zero, CEP/1.1774) 

Hieght cf burst and yield are selected from univariate 
normal distributions. 

HOB = N (desired HOB, PEH/.67) 

Yield = N (nominal yield, nominal vield/10) 



F. ASSESSING DA23AGE 

Fcr damage assessment it is assumed that HOBcpt was 
selected as though overpressure was the governing effect. 
This is dene to avoid reading Figs. 2.1 and 2.2 into the 
data base. This is dene for two reasons. 

1. The data in the figures is very coarse and should not 
be accepted as highly accurate [Ref. 1: preface], 

2. A linear function results from the assumption. 
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If the assumption that the weapon will function at HOBcpt 
for overpressure is made, then one variable, the HOB, can be 
omitted. If the employed weapon functions at HOBcpt then 
the scaled HOB is simply the HOBopt for a 1 KT weapon. 

Using this HOB and the ranges found in Fig. 2. 1 for 
desired c ver pre ssures and in Fig. 2.2 for desired dynamic 
pressures, it is possible to plot yield against range with 
pressure as a parameter. Figs. 4.5 and 4.6 show these 
results. 



1 




Figure 4.5 Overpressure from Optimum HOB. 



If the weapon yield and the required overpressure or 
dynamic pressure for a particular target defeat is known, it 
is a simple mar ter for the computer to calculate the range 
at which defeat will occur. An example will demonstrate 
this. 
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Figure 4.6 Dynamic Pressure from Optimum HOB. 




Given: A de-conation results in a 10 KT yield. An AFC is 

damaged cy 15 psi overpressure. 

Find: The distance from ground zero at which the AFC is 

damaged. 

Solution: From Fig. 4.5 locare 10 KT on rhe lower scale. 

Go up to the 15 psi line. Sead 790 meters from the left 
scale . 



In rhe model. Eg. 4.1 is used to determine the range. 
The coefficient A is a function of overpressure or dynamic 
pressure and is determined from Figs. 4.5 and 4.6 when the 
yield is 1 KT. 



RD = A* Yield 17 “ 



(4.1) 
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The same simplification can be made concerning arrival 
time. Equation 4.2 is the result. 

T = {(.0028-R) - . 45} *Yield 1/3 ( 4 -2) 

It must be pointed out that this method will predict a 
larger HD fcr overpressure than is correct when overpressure 
is not the governing effect. However, overpressure is 
almost always the governing effect. (Radiation is usually 
the governing effect only against personnel in the open.) 
But, when radiation is the governing effect, enhancing the 
effects of overpressure to their optimum will still not 
outweigh the effects of radiation. A soldier killed by 
radiation will still be killed by radiation. The only error 
is that a soldier who survives radiation and should have 
survived overpressure will now have a decreased chance of 
survival. 

The alternative to the simplified method is to revert to 
Eq. 2.1 with a table look-up and 2 way interpolation. 

G. VARIABILITY OF EFFECTS 

An object positioned at a certain range from a given 
yield will be subjected to overpressure, dynamic pressure, 
thermal radiation and nuclear radiation. In many cases one 
damage mechanism will dominate the others. It would be a 
simple matter to use only the dominant mechanism in 
determining damage. In an expected value model this 
technique is entirely justified but such a proceedure is not 
followed in this model. 

In a probabi listic model, such as this one it is 
entirely possible for a target element to "beat the odds" on 
cverpressure and be killed by secondary missiles brought 
about by dynamic pressure. If there was a 75& chance of 
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damage by overpressure and only a 103 chance of damage by 
dynamic pressure in a 503 threshold model, the element would 
be killed by overpressure every time. In this model, each 
mechanism gets a Monte Carle chance to cause damage. In the 
above example the subject has a 77.53 chance of damage. 
More importantly, it has a non-zero chance of survival. 

If the ED from Fig. 4.5 is found to be 700 meters (10 
psi and 10 KT) then what are the probabilities for survival 
at a closer range and for damage an a farther range? The 
expected value model would set these probabilities at zero. 
This model will allow survival at a closer range and damage 
at a farther range. To do this requires a probability 
distribution. 




Figure 4.7 Variability of Effects. 

Fig. 4.7 [Ref. 4: p. D-3] is used in the FM 101-31 
series and is followed in this model. The 203 variance 

curve is the one which is used. Plotted on normal plotting 
paper, it locks like Fig. 4.3. 
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Figure 4.8 Standard Deviation is 20% of R/HD. 
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Kill 



According to Fig. 4.7, a -target element at ground zero 
has some finite chance of survival while an element at an 
infinite range may still perish. This is a logical result, 
of the normal probability assumption. To approximate 
reality, any element with less than 0.5% chance of survival 
is killed while anything with more than 99.5% chance of 
survival always survives. Thus, only when R/RD is between 
.5 and 1.5 is the Monte Carlo chance taken. An example will 
demonstrate how the method works. 

GIVES: A 6 KT burst. Tanks are damaged by 20 psi 

ever pressure . 

FIND: Damage to tanks at 200 and 900 meters from GZ. 

SOLUTION : From Fig. 4.5 the range at which a tank is 

damaged is 500 meters. Correctly stared, this is RD, the 

distance at which the tank has a 50% chance of survival The 
first tank has a distance/RD of .40. From Fig. 4.8, the 
probability of survival is less than 0.5% and the tank is 
killed . 

The ssccnd tank has a distance/RD of .555. A uniform 

(0,1) draw is made for the tank and the number of .60 is 

drawn. Since .555 is less than .60, the tank survives. Had 
.554 cr less been drawn, the tank would have been killed. 

H. TARGET RESPONSES 

Target elements are of two types - tanks and personnel. 
However, any target element which is damaged by blast, 
thermal radiation or nuclear radiation may be easily input 
into the model. Such elements would include APC's, trucks, 
buildings, aircraft, forests, bridges, etc. 

In this model the tanks are either fully functional or 
killed. Personnel may have an additional state - dying. 
This is caused by the absorption of nuclear radiation, 
while the various mechanisms may kill a soldier upon burst. 
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only the nuclear radiation will cause a latent, lingering 
death. Fig. 4.9 [Bef. 4: p. C-9] is used to classify 

radiation casualties into the states of: 

1. Dead. 

2. Incapacitated and dying. 

3. Impaired and dying. 

4. Safe and fully functional. 



PHYSICALLY DEMANDING 




Figure 4.9 Personnel Besponses to Radiation. 

The model simplifies Fig. 4.9 into Fig. 4.10 with the 

following : 

1. Anyone receiving 18,000 rads is declared dead 
imme diately . 

2. Anyone receiving 8,000 to 18,000 rads is immediately 
and permanently incapacitated. He cannot shoot, move or 
communicate but is still alive and should be considered 
a potential target. 

3. Anyone receiving 3,000 to 8,000 rads is immediately 
incapacitated but in time is upgraded to functionally 
impaired. The time is uniformly distributed from 30 to 
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45 minutes. Eventually the person dies. Death occurs 
in a uniformly distributed amount of time from 2 to 5 
days. 

4. Anyone receiving 2,000 to 3,000 rads is immediately 
placed in the functionally impaired state where he 
remains until death occurs in 2 to 5 days. 

5. Anyone receiving a dose greater than his own personal 
lethal limit but less than 2,000 rads is unaffected 
until some latent period has passed at which time he 
becomes functionally impaired. This time is uniformly 
distributed from 30 to 60 minutes. Eventually, the 
person dies in 2 to 5 days. The lethal limit is 
normally distributed about 650 rads with a standard 
deviation of 50 rads. 

6. Anyone receiving less than his own personal lethal 
limit is unaffected. 

I. OUTPUT 

When the model is conpletely finished it will print all 
of the pertinent information concerning each target and 
employed round. Final output will be; 

1. Burst parameters of GZ, HOB and yield for destruction 
of each target. 

2. Location and disposition of each element in each 
target. For personnel not immediately killed, the 
incapacitation and/or impairment times will be listed as 
appropr iat e . 
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Figure 4.10 Simplified Personnel Responses to Radiation. 



J. SUMMARY 

This model will, given several nuclear capable artillery 
batteries, friendly locations and enemy locations, select 
optimum choices of batteries ana yields to fire at each 
enemy target. It will, when necessary, allow for troop 
safety fcy offsetting the aim point of the weapon. It will 
simulate delivery of each round and stochastically assess 
damage to each individual element within the target. Tanks 
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will t s declared either alive or dead while personnel will 
fee allowed a dynamic transition from living to dead due to 
nuclear radiation. 
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V. MODEL SOUTINES 

This chapter will explain the important sections of the 
actual model in detail. Actual lines of the SIMSCRIPT model 
(printed in capital letters) will be presented followed by 
an explanation of what they do or mean. The complete modal 
is included as Appendix B. Each routine and function is 
discussed as a sepatate section in this chapter. 

In reading about each routine and function it may be 
helpful tc refer to Fig. 5.1 to keep the flow of the model 
in order. 

A. P8EAMELE 

The preamble extabiishes global variables, identifies 
permanent and temporary entities and their attributes, 
groups entities into sets and identifies user defined 
funct ions . 



PERMANENT ENTITIES 

EVERY EATTERY HAS A SIZE, AN XB, A YB AND SOME NOC. ROUNDS. 
EVERY COMPANY HAS A XC , A YC, A ZC AND A CO. RADIUS 
EVERY TARGET HAS A TYPE, A XT, A YT, A ZT, A RT 

AND MAY OWN A TANK. SET, A TROOP. SET, AND A LISTING 

The EATTERY is an entity representing a nuclear capable 
artillery battery. The SIZE is either short range cannon 
(155mm, SIZE = 1) or medium range cannon (8 in., SIZE = 2). 
XB and YE are grid coordinates. NUC. ROUNDS is the quantity 

cf nuclear rounds the BATTERY has available to fire. 
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Figure 5.1 Routine/Function Routing Chart. 

The COMPANY is an entity representing a friendly unit. 
It may be an individual soldier or vehicle cr a collection 
cf such. Its reason for existing is to cause offset 
targeting due to friendly safety constraints. XC, YC, ZC 
are the grid coordinates and altitude. CO. RADIOS is the 
radius of the unit. 
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The TARGET is an entity which is to be attacked. XT, YT 
and ZT are grid coordinates and altitude. RT is the target 
radius cr equivalent as explained in Chapter III. The 
TARGET owns three sets. They are sets of TANKS, TROOPS and 
a LISTING of feasible SOLUTIONS. 

TEMPORARY ENTITIES 

EVERY TANK HAS A X . TANK, A Y . TANK, A TANK. DEAD, 

AND MAY BELONG TO A TANK. SET 
EVERY TROCP HAS A X. TROOP, A Y. TROOP, A TRANS . FACTOR, 

A TROOP. DEAD, A LETHAL. DOSE, AN ACCOM. DOSE, 

A T.LETH, A T. IMPAIR. NUC 
AND MAY BELONG TO A TROOP. SET 
EVERY SOLUTION HAS A U NIT . TO . FIRS, AN XDGZ, A YDGZ , 

AN IYIELD, A JROW, A KCCL AND A PCT 
AND BELONGS TO A LISTING 

The TANK has coordinates X.TANK and Y.TANK. The 

altitude is assumed equal to that of the target to which it 
will te assigned. TANK. DEAD is an integer number, either 0 
cr 3 . AO means the tank is fully functional. A 3 means it 
has suffered a catastrophic kill. 

The TRCOP has coordinates X. TROOP and Y. TROOP. The 

altitude is assumed equal to that of the target to which it 
will be assigned. The TRANS. FACTOR is the radiation 
transmission factor and is currently either a 1 for exposed, 
.3 for protection inside an APC and .001 for a fcx hole. 
TROOP. DEAD is similar to TANK. DEAD except for two additional 
values. A 1 is functionally impaired and a 2 is 

incapacitated. LETHAL. DOSE is that radiation dose above 
which the TROOP will eventually die. ACCUM.DOSE is the 
radiation dose the TFCOP has currently accumulated. T.LETH 
is the time when TROOP. DEAD is to be set equal to 3. 
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T. I M P AI B . NUC is the time when TROOP. DEAD is to be set equal 

to 1 . 

A SOLUTION is an answer to the nuclear target analysis 
problem. Its attributes are the firing unit, desired ground 
zero, expected coverage and indicies to locate the table, 
row and cclumn where the solution was found. It is filed in 
a LISTING belonging to a TARGET. The LISTING is ranked 
first by high coverage and then by low yield. 

GLOBAL VARIABLES 

There are two sets of equivalent arrays. One pertains 
to the SRC and the ether to the MRC . The SRC has 9 rows 
while the MRC has 19 rows. In the following, only the SRC 
is explained but the description extends equally to the SRC. 

SRC. P .COVER AGE is the 4 dimensional matrix which 
contains the coverage probabilities. The first subscript is 
the yield. The second is the type of target. The third and 
fourth are the row and column of the table described by 
yield and type. The current model allows seven types of 
targets. They are; 

1. Exposed personnel. Immediate permanent casualties. 

2. Exposed personnel. Immediate transient casualties. 

3. Personnel in fox holes. Immediate permanent casualties. 

4. Personnel in fox holes. Immediate transient casualties. 

5. Moderate damage to Tanks. 

6. Moderate damage to Wheeled Vehicles. 

7. Moderate damage to Artillery. 

SRC. RANGES is the array of ranges that the system is 
capable of. (See Fig. 3.5) 

SRC. TABLE. RT is the 3-d imensicnal matrix which contains 
the RT fer the various tables. The first subscript is the 
yield. The second is the type of target. The third is the 

column of the table. (See Fig 3.5) 
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SBC. YIELDS is the array of possible yields. Each weapon 
has two yields available. The SRC has a . 2 KT and a 1KT 
while the MEC has a 2 KT and an 8 KT. 

MSB. SEC is a 2- dimensional marrix which contains values 
of minimum safe distances. The first subscript is the row 
of the table and the second is the yield. MSD is only a 
function cf weapon size, range and yield but nor target type 
or degree of damage. 

CD90.SRC, CEP. SEC and PEH.SRC are arrays with only rhe 
range as a variable. They are a function cf rhe delivery 
system only. The subscript is the row of any table for rhe 
applicacle weapon system. 
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variables and call subroutines and funcrions which do the 
bulk cf the work. 
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POLICY is an array which represents the commander ' s 
guidance cn each of the types of targets. Its length is 
NTABLES * one for each of the seven types of targets. If 
the predicted coverage on a target does not equal or exceed 
the POLICY for that type of target, then the solution is 
abandoned. 

RANGE is the range from the battery to the target for 
which a solution is being sought. 

The remainder cf MAIN is involved with calling 
subroutines to solve the problem. MAIN first calls 
HEAD. RESERVE which reserves and reads all of the global 
variables. MAIN then calls INITIALIZE. It then loops 

through every TARGET and every 3ATTERY in an attempt to find 
adequate SOLUTIONS . To do this it calls ANALYSIS which, in 
turn, calls other subprograms. After returning from 
ANALYSIS, the MAIN is ready to employ the weapons. It 
removes the first SOLUTION from the LISTING of each TARGET 
and calls DETONATION which assesses damage. Upon returning 
from DETONATION , the MAIN lists the attributes of the 
permanent entities and terminates. 

C. SUBROUTINE READ. RESERVE 

This subroutine reserves and reads the global variables 
as defined in the PREAMBLE 

D. SUBROUTINE INITIALIZE 

This is a surrogate for a driver program. As stated in 
the abstract, this program is intended as a subroutine for a 
larger simulation - one which already has entities such as 
TARGET, TANK and TROOP. Since these entities do not exist, 
INITIALIZE creates them. A detailed description of the 
subroutine will not be offered hers since it is unimportant 
and easy to read. It will suffice to say that it uniformly 
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distributes TANK and TROOP entities within the RT of the 
TARGET while assigning attributes to the temporary entities. 

E. SUBROUTINE ANALYSIS 

This subroutine is the origin of SOLUTION. It is here 
that a decision is made to create a SOLUTION and then to 
file it in a LISTING or discard it. The first action is to 
determine the size cf the EATTERY which has been passed to 
it. It then determines if the TARGET is within range. If 
not, it returns to MAIN without a SOLUTION. If the range is 
less than the maximum weapon range, then it finds the proper 
row (J) to use as an entry argument. It then loops ever all 
yields to determine coverage for each yield. For each yield 
it determines the proper column (K) to use as the ether 
entry argument. If it cannot find any (RT too big) it loops 
to the next yield cr returns to MAIN if if was at the 
largest yield. 

If a row and column were found, then a SOLUTION is 
created and the index method is used to compute the initial 
coverage. If the ST is less than or equal to the first 
listed in the table, then the coverage in the first column 
is reported. If not, the function EXP. COVERAGE is called. 
This returns the expected initial coverage using the index 
method. 

A decision is made concerning the coverage and the 
command guidance (POLICY) for that type of target. If the 
coverage is adequate, then OFFSET. COVERAGE is called. This 
will return with a new coverage which is less than cr equal 
to the cne returned by EXP .COVERAGE. If this new coverage 
is still adequate, the attributes of SOLUTION are determined 
and it is filed in the LISTING cf the TARGET for which it 
was called. Control is then returned to MAIN. 
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F. 



SOBBOOTINE EXP. COVERAGE 



This function was originally established to perform the 
interpolation with the index method. However, as it became 
necessary tc perform interpolation in other subprograms, a 
dedicated interpolation function (INTER?) was written. 
Currently this function simply determines the size of the 
BATTERY and properly formats the arguments for a call to 
INTER?. It then returns with the expected coverage to 
ANALYSIS. 

G. SOBBOOTINE OFFSET .COVER AGE 

This subroutine is called when the initial analysis is 
at least adequate. It determines if any offset is necessary 
and, if so, performs an analysis using the numerical method. 

The first order cf business is to call MAX. DISTANCE 
which actually locates the new desired ground zerc. If 
there is no displacement, then there is no change in 
coverage and control is returned to ANALYSIS. MAX. DISTANCE 
will also return a flag which may mean nc feasible solution 
exists and coverage is set equal to zero. This will cause 
the solution not to he filed in a LISTING. 

Next, the battery size is determined and the entry 
arguments cf Fig. 3.7 are calculated. Depending on the 
value of d/CD90 (Y2) ere cf three things can happen; 

1. Consider the target as a point target. 

2. Analyze as instructed. 

3. Degradation is insignificant. Return. 

Currently, a point target analytical technique is not 
included, tut the technique is slightly easier than this 
one. The only time an area target would be analyzed as a 
point target is if the RD is at least 19 times greater than 
RT. This only happens with large yields. Case 2 is the 



required (assuming 



only one for which further analysis is 
case 1 dees not happen) . 

Case 2. d/CD90 is located on the ordinate and (with a 
call tc INTERP) the measured distance is found. As 
explained in Chapter III, the measured distance (DELX ) is 
multiplied by the original CD90/RT and used to find the 
predicted coverage due to offset. A double inter polat icn is 
necessary to find the coverage. Control is then returned to 
ANALYSIS. If the coverage is still adequate, the SOLUTION 
is filed in a LISTING. 

H. FUNCTION INTERP 

This is an elementary interpolation function, is self 
explanatory and is net further explained here. 

I. SUBROUTINE MAX. DISTANCE 

This subroutine performs the approximation to the 
non-linear programming problem as stated in Chapter III. 
The lcgic fellows the graphical solutions as shown in Figs. 
4. 1 and 4 . 2 

J. SUBROUTINE DETONATION 

The purpose of this subroutine is to assess damage new 
that the analysis is complete. The best SOLUTION for a 
TARGET has teen removed from its LISTING and passed to this 
subroutine. 

Again, a decision is made regarding the BATTERY sire. 
The standard deviation for ground zero location is found Eq 
3.1. The PEH and desired HOB are selected from the tables. 

Tc determine ground zero the radial miss distance, R, is 
computed as N(0 , CEP/ 1 . 17 74) while an angle. A, is picked 
from U(0,2 tt ) . Ground zero is then placed at the polar 
coordinates of (R,A) 
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The actual yield is selected from a normal distribution 
with nominal yield as the mean. (Actual variances are 
classified) The standard deviation is assumed to be cne 
tenth cf the mean. 

HOB is selected from a normal distribution with DHCB as 
the mean and PEH/.67 as the standard deviation. 

The actual GZ, yield and HOB have now been fixed and the 
burst parameters are output by the statement; 

LIST X, Y, YIELD, HOE 

EACH TANK in T ANK . SET (TARG E T) and TROOP in 

TROOP . SET (TARGET) is examined to check for damage. As 
stated, the TANK can only be fully operatinal (TANK. DEAD = 
0) or killed (TANK. DEAD =3). The only mechanism which has 
any reasonable chance to destroy a tank is overpressure. 
The overpressure at which a TANK has a 503 chance of 
survival is 25 psi. Thus; 

LET RD. OVER. PRESSURE = OVER . PRESSURE (YIELD, 25) 

will return the distance from GZ where 25 psi will be 
encountered and will be used to determine if the TANK is 
destroyed . 

IF (R/RE. OVER. PRESSURE LE .5) OR 

(.01* (NSCTRCN(R, YIELD) PLUS GAMMA(R, YIELD) ) GT 18000) 

is a compound IF statement which destroys the TANK 
(TANK. DEAD = 3) if either of the following conditions are 
met . 

1. If the probability of destruction by overpressure is 
greater than 99.5% or 

2. If the crew inside would receive more than 18,000 
rads of radiation. 
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The TRCCP is put through all four damage mechanisms 
before being declared a survivor. The thermal level at 
which a TROOP has a 50% chance of survival is 20 cal/cm . 
Thus ; 

LET RD.HEAT = HE AT ( yield ,20) 

will return the distance from GZ where 20 cal/cm will be 

encountered and will be used to determine if the TROOP is 
killed. Tc become a thermal casualty a compound IF 
statement must be true: 

The TROOP must be exposed (TRANS . F ACTOR (TROOP) = 1) and 
The TRCCP must receive enough thermal energy to become a 
therial casualty (PROB ABILITY (R/RD .HEAT) LT CHANCE) and 
His reaction time must be greater than . 1 second. 

If the TROOP survives the first test, he may become a 
casualty due to overpressure. The overpressure at which a 
TROOP has a 50% chance of survival is 10 psi. Thus; 

LET RC. OVER. PRESSURE = O VER . PR ESS URE ( YIELD , 1 0) 

will return the distance from GZ where 10 psi will be 
encountered 

and will be used to determine if the TROOP is destroyed. 

If R/RD is less than .5, the TROOP is killed (TROOP. DEAD = 
3). Otherwise he must take a Monte Carlo chance; 

IF PROBABILITY ( R/RD . OVE R. P RESSURE) LT CHANCE 

If the above statement is true then the TROOP is killed. 

If the TROOP is still alive, he must now survive the 
effects cf dynamic pressure. This section is identical tc 
overpressure with the exception of a call to 
DYNAMIC. PRESSURE instead of OV ER. PRESS URE to get 

RD. DYNAMIC. PRESSURE 
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After all of this, if the TROOP is still alive, nuclear 
radiation gets its chance. The combined effects of neutron 
and gamma gamma radiation are added to the current radiation 
level . 

ADD TRANS. FACTOR(TRCOP) *(( NEUTRON (R , YIELD) +GAMMA (R, YIELD)) ) 
TO ACCUM. DOSE(TROOP) 

If the TROOP has accumulated more than 18,000 rads he is 
killed (TECCP. DEAD = 3). 

If the ACCUM. DOSE is between 8,000 and 18,000 rads he is 
placed in the incapacitated stats <T ROOP. DEA D (TROOP) = 2) 

until he is killed in 2 to 5 days. 

LET T. LETH (TROOP) = R ANDI .F (2*24*60, 5*24*60, 1) 

If the ACCUM. DOSE is between 3,000 and 8,000 rads he is 
placed in the incapacitated state fcr 30 to u5 minutes. 

LET T. IMPAIR. 1JUC (TROOP) = H A NDI. F (30 , 45 , 1 ) 

At the end cf this time he is placed in the functionally 

impaired state (TROOP. DE AD (TROOP) =1). He is killed in 2 
to 5 days. 

If the ACCUM. DOSE is between 2,000 and 3,000 rads he is 
placed in the functionally impaired state until he is killed 
in 2 to 5 days. 

If the ACCUM. DOSS is greater than what he can tolerate 

(LETHAL. DCSE (TROOP) ) but less than 2,000 rads the TROOP is 

unaffected for a short period of time - 30 tc 60 minutes. 

LET T. IMPAIR. NUC (TROOP) = R ANDI. F (30 , 60 , 1) 

At the end cf this time he is placed in the functionally 

impaired state. He is killed in 2 to 5 days. 

If ACCUM. DOSE is less than LETHAL .DOSE (TROOP) then he is 
fully functional and will not be killed as a result cf the 
current detonation. 
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K. FUNCTION 07EB. PHESSUBE 

This function determines the distance from GZ at which a 
target element will encounter a 50% chance of survival for 
the achieved yield and required oper pressure. 



L. FUNCTION DYNAMIC. EHESSURE 

This function determines the distance from GZ at which a 
target element will encounter a 50% chance of survival for 
the achieved yield and required dynamic pressure. 



B. FUNCTION ABRI7AL.TIME 

This function executes the arrival time function and 
returns the time after detonation when the blast wave will 
pass the given distance fromn GZ. 

N. FUNCTION GAMMA 

Equation 2.4 is executed by this function. For the 
achieved yield and distance from GZ the amount of neutron 
radiation is returned. 



0. FUNCTION NEUT20N 

Equation 2.3 is executed by this function. For the 
achieved yield and distance from GZ the amount of neutron 
radiation is returned. 



P. FUNCTION HEAT 

Equation 2.2 is executed by this function. For the 
achieved yield and thermal requirement the distance from GZ 
where a 50% chance of lethal burns occurs is returned. 
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Q. FUNCTION PROBABILITY 



This function approximates zhe area under a N<0,1) 
probatility curve. It is a table lock-up for Fig. 4.8. 
Only entry arguments (R/RD) of .5 and greater are passed to 
it. If the argument is greater than 1.5 zhe probabilitiy is 
set tc 0. 
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VI. ENH ANCEM ENTS AND EXTE NSIONS 



A. OPTIMOM ASSIGNMENT OF BATTERIES. 

This model assigns batteries to fire on an essentially 
random basis. There is no thought given to optimizing the 
coverage. Choosing the highest coverage for each target or. 
a first in first out basis may not provide maximum total 
coverage. An assignment problem must be solved to do this 
correctly. 

Table IV is a potential nuclear load for three batteries 
in the following example. 



TABLE IV 

Nuclear Load by Battery 



Battery Type 



1 

2 

3 



ssc 

SRC 

MRC 



No . cf rounds 
,2 KT 1 KT 2 KT 8 KT 



1 

1 

0 



2 

1 

0 



0 

0 

1 



0 

0 

2 



Table V shows the coverages that can be obtained with 
the various batter y/yiel 3 combinations or. five targets. A 
zero means that the target is either out of range or 
coverage is less than command guidance. Note that battery 1 
has 2 rounds of 1 KT yield and, thus, has 2 identical rows 
for the 1 KT rounds. A similar situation is true for the 8 
KT rounds of battery 3. 
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TABLE V 

Target Coverage 



Ta rge t 

Battery Yield 12345 





.2 


KT 


.72 


.00 


.00 


.00 


.00 


1 


1 


KT 


.82 


.80 


.00 


.45 


.00 




1 


KT 


.82 


.80 


.00 


.45 


.00 




. 2 


KT 


.75 


.00 


.42 


.00 


.00 


2 


1 


KT 


.65 


.00 


.51 


.42 


.00 




2 


KT 


.92 


.63 


.47 


. 97 


.92 


*4 


8 


KT 


.57 


.00 


.00 


.00 


.00 




8 


KT 


.57 


.00 


.00 


.00 


.00 



Tc assign battery/yield combinations to 



tar cet 



according 

following 


to the 

• 

f 


current te 


Target 


Eannery/Yi 


.elc 


Covera 


1 


3/8 


Kt 


. 97 


O 


1/1 


Kt 


. 30 


3 


2/1 


Kt 


.51 


4 


3/2 


Kt 


. 97 



This arangement results in a tonal coverage of 3.25, cun 
target number 5 cannot be engaged because the only feasible 
round for it was employed on target number 4. The 
assignment problem that needs to be solved is one that will 
maximize total coverage. This case is very easy tc solve by 
inspection . 
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The cpti 


mum assignment 


is; 


Target 


Hatter y/Y ield 


C overa 


1 


3/8 Kt 


. 97 


2 


1/1 Kt 


.80 


3 


2/1 Kt 


. 51 


4 


1/1 Kt 


.45 


5 


3/2 Kt 


.92 



The optimum total coverage is then 3.65. 

In the current model, only the percentage of coverage is 
used to elvaluata the value of the target. This may net be 
the item to maximize if the various targets are not of the 
same type and population. For example, a 97% coverage on an 
infantry company is not as valuable as a 50% coverage on a 
tank regiment. To make the matrix in Table V mere 
realistic, a variable such as; 

Value = coverage * element value * target population 
should be used. 



B. VARIATION OF MS D 

Typically, an operations order specifies the troop 
safety guidance as a negligible risk to unwarned exposed 
personnel. This is the category used in the model tc 
determine MSD. Fig. 6.1 shows the different separation 
distances that must be achieved for the safety categories 
C- , 

1. Unwarned exposed. 

2. Jiarned exposed. 

3. Earned protected. 
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Each category has the 3 levels of radiation exposure 
status (RES) of ; 

1. Negligible. 

2. Moderate. 



3. Emergency. 







UNWARNED 




WARNED 




WARNED 








EXPOSED 


EXPOSED 


PROTECTED 




RANGE 


NEC 


MOD 


EMER 


NEG 


MOD 


EMER 


NEG 


MOD 


EMER 






7500 


6990 


*♦500 


44QU 


4100 


?2Jr 


3400 


3000 


1600 




JOO 0 


750 0 


0900 


*♦500 


44J J 


4100 


33 JO 


3500 


3 100 


160(1 




4U0O 


7500 


7000 


*♦5^0 


44 Ou 


4100 


3300 


3500 


3100 


1700 




‘jGljO 


7bC0 


70u3 


46 JG 


4500 


4200 


3300 


3600 


3P0C 


1700 




o0 jO 


7b00 


7C00 


460C 


4500 


4P90 


3300 


3600 


3P00 


1 70 r 




7u j0 


/6C0 


7U00 


46JU 


450* 


4 POO 


34 Ju 


3700 


33or 


1000 




3JjO 


7oPU 


7100 


*♦60 U 


45 Jo 


4200 


34 0( 


3700 


33f 0 


1P00 




9000 


77C0 


71u0 


*♦700 


4bOU 


4 300 


34UC 


3500 


3400 


1 POO 




IOOuO 


/700 


7100 


*♦700 


460b 


4300 


3400 


3^00 


3500 


1P00 




IIQ'JO 


7700 


7 luO 


*♦700 


460U 


4 300 


3500 


4000 


3*oe 


l«0O 




IcOoO 


7700 


7<i00 


4 70u 


4600 


4300 


350f 


4000 


3600 


1°00 




i jujQ 


7bOO 


72uQ 


'♦aoo 


470U 


4400 


35UC 


4100 


370 0 


1°00 




*4G W 0 


7aoc 


7200 


9BJU 


4700 


4400 


3500 


4 2c 0 


3700 


1°90 




1.50^0 


7b U0 


7*00 


*♦390 


470o 


4400 


36 JC 


4300 


3 of 0 


2000 




) out 0 


7*00 


7300 


*♦000 


470 U 


4400 


36 JC 


4300 


3oou 


PO00 




*. 70 jO 


7900 


7300 


490b 


4800 


4500 


36 J(J 


4400 


3°O0 


200C 




1 d0 o 0 


7^00 


73o0 


*900 


4H00 


4500 


36J( 


4600 


4000 


2000 




Figure 6 


.1 


Minim um 


Separation Di 


stances. 


In Fig. 


6.1 


th 


S MSD 


for a 


. negl 


igible 


ri sk 


to unwarned 


exposed 


and warned exposed 


personnel 


ar 


e 7700 and 4600 


meters r 


aspect 


ive 


iy. 


Recall 


. that 


MSD is 


the reason for 


offset 


and o 


f fs 


et i 


— 

O u 


he reason 


fo 


r some 


very severe 



degradation in coverage. If the 4600 meter MS D could be 



used, then the displacement could be reduced by 
This could make a significant difference in 
outcome of the simulation. 

Expansion of the MSD category would greatly 
use of the model as an analytical tool. 



430C meters, 
the final 

expand the 
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C. COLLATERAL DAHAGE AVOIDANCE. 

Chapter 14 of Ref. 2 contains least separation 
dist ances (LSD) to avoid damage to many entities such as 
buildings, forests, bridges and civilian populations. This 
imposes the same type of offset requirement on the aimpoint 
as USD dees except that it precludes damage to entities 
ether than friendly troops. 

Use of this idea would require implementation cf a 
routine similar to that which determines offset except that 
the description of what is being avoided must also be 
included . 

D. HUITIELE TARGETS 

As it has already been shown, the center of a target 
need not be the aimpeint. It would be highly advantageous 
to destroy more than cne target with a single burst. Fig. 
6.2 shews 3 targets and their associated maximum 



o 



Target 1 



Figure 6. 2 




Multiple Targets. 



displacement distance to insure minimum acceptable coverage. 
The "triangular" area that is inside all circles is a 
feasible area for each target. That is, a round employed 
anywhere in the "triangle" will destroy all 3 targets. 
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A method to evaluate groups of targets would greatly 
expand the use of the model as an analytical tool. 

E. F ALLCOT 

Generally, nuclear weapons are detonated at or above an 
altitude that provides a 99% assurance of no fallout (Eqs. 
3.4 and 3.5). However, fallout may be an effective barrier 
to movement or cause mass evacuation of a contaminated area. 
Control of fallout after employment is impossible due to 
atmospheric conditions and the user may become the victim. 

Should a ground hurst occur, it is necessary to predict 
and simulate fallout. Predicting fallout requires that 
meterclogical data at altitudes up to and including cloud 
top must be known. A deterministic model for predicting 
fallout is relatively simple but actually simulating fallcut 
may require assistance from experts in meteorology. 



F. ACTIONS TAKEN DOBING AND AFTER FALLOUT. 

While the detonation of nuclear weapons will certainly 
cause seme unique actions and reactions, there are some 
specific actions required of soldiers upon encountering 
fallcut. 



1. C_ptimum Time cf Exit. 



For perscnnel caught in a fallout area (presumably 
in covered fox hcles cr shelters) there is an optimum time 
to exit the area to minimize the absorption of nuclear 
radiation. 

While soldiers stay in the protective confines of 
their fex hcles, they take advantage of the shielding. They 
also prolong their exposure time. To immediately leave 
would reduce their exposure time but would expose them tc a 
much higher dose rate outside their shelters. 
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For given initial dose rate, decay rate, shielding 
and exposure time to a safe area, there is an optimum time 
to leave. 

2. Crossing Prob lem . 

Given that a unit must cross a radiation area, there 
is a not-eariier-than time before which they may not enter. 
It is a function of 5ES , fallout pattern and decay rata, 
speed of crossing and shielding during the crossing. 
Inclusion of this type of problem into the the maneuver 
section of a driver routine would certainly have an impact 
upon the outcome of the simulation. 

3* Deco nt a m inat ion . 

Actions taken to decontaminate vehicles and 
personnel are relatively simple if only transportation to 
and from a decontamination site and time taken to 
decontaminate are considered. The main problem with 
decontamination is one of supply. 

G. ADDITIONAL TARGET ELEHE8TS. 



The 


current 


mod 


el has 


only 


tanks and troops i: 


i the 


target areas. 


Any 


target 


may 


fce 


added to the target 


array 


by defin: 


Lng it 


as a 


tempor 


ary 


err 


tity and allowing a target 



to own it. To damage the entity requires the modeler to 
call one or more of the damage functions with the 
appropriate arguments and then decide if the entity is 
damaged. The procedure is virtually identical to what is 
currently used. 



84 



APPENDIX A 



BEGBESS ION ANALYSIS 



A. GAMMA B ADI AT ION 

Since the model needs the gamma radiation level at a 
specific range for a given yield, it is advantagcus to 
raplct Fig. 2.7 with dose as the dependent variable. Fig. 
A. 1 is this new plot. 




Surrendering this plot to the APL linear regression 
package for the equation; 

In (Dose) = A + 3 * Range + C * In (Yield) 

fails tc provide an adequate solution. Therefore, the model 
was solved in pieces and reassembled. 
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The procedure was to regress each line independently as 
In (Dose) = A + B * Range 

for each yield. The coefficients A and B are obtained for 
each yield and are listed in Table VI. 



TABLE VI 

Coefficients A and 3. 



Yield A B 



1 


11.18 


-.0060 


2 


11.56 


-.0056 


3 


11.76 


-.0054 


4 


11 .91 


-.0052 


5 


12.02 


-.0052 


5 


12 . 22 


-.0051 


7 


12.22 


-. 0050 


8 


12.33 


-.0050 


9 


12.38 


-.0049 


10 


12.42 


-.0049 



Fcr each yield, the coefficients are plotted as shewn in 
Figs. A. 2 and A. 3. 




86 



-.005 



□ 
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-.006.- 
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Figure A. 3 Begression Coefficient B (Data). 



The coefficients can now be regressed against yield. From 
Fig. A. 2 the first coefficient appears to vary as the square 
root of yield. Therefore, the equation; 

A = C 0 + C 1 -/Yield 
is fit with the following results. 

C 0 = 10.74 

C = .5599 

This equation is plotted in Fig. A. 4 along with the 
original coefficient data. The fit is qurte good. 
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There is a temptation to say that the coefficient, E, is 
really a constant since the curve is very flat and the 
represented values are so small. However, it must be 

remembered than this is to be multiplied by the range and 
the result used as an exponent. Any small error may have 
catastrophic results. From Fig. A. 3 a straight line of the 
form ; 

E = C 0 + Cj * Yield 
is fit with the following results. 

C 0 = ” .0053 

C = .0001 

i 

The equation is plotted in Fig. A. 5 along with the original 
coefficient data. 




Assembling the equation results in the following 
equation . 



G 



46166 



.5599/Yield . 0001 -Yield *R 

e e 

. 00 5 8 • R 
e 



(A.1) 



68 



A. 6 for even yields 



Equation A. 1 is plotted on Fig. 
along with the original data. 




It is important that the equation be accurate in the 
vicinity cf 650 rads as this is the mean of the lethal dose 
distribution and is the dividing point between surviving and 
dying. 

E. NEUTECN RADIATIOS 

Like gamma radiation, the modal needs the neutron 

radiation level at a specific range for a given yield. 
However, Fig. 2.6 appears to have straight lines which can 
be exploited in a regression analysis. Fig. 2.6 has been 
converted to a metric scale and is presented in Fig. A. 7. 

The neutron regression model was also split into two 
submodels. Each line was regressed independently as; 

Ranee = A + B * In (Yield) 



for each dose. 



39 




Yield (KT) 

Figure A. 7 Heutron Radiation. 

The coefficients A and B are obtained for each dose and are 
listed in Table VII. 
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TABLE VII 
Coe f ficients 



1 



Dose 


A 


B 


30 


1272 


171 


100 


1093 


158 


300 


9 18 


156 


1000 


766 


148 


30 CO 


641 


129 


10000 


4 61 


138 




F cr each dose, the coefficients are plotted as shewn in 
Figs. A. 8 and A. 9 . The coefficients can now be regressed 
against the dose. 




Figure A. 8 Begression Coefficient A (Data). 



For the first coefficient. A, the nodal is 



with 



In (A) = C 0 ♦ Ci * In (Dose) 



C = 7.77 

o 

Ci = -.16 96 



9 1 



180 - 




Fcr the secor.d coefficient, B, the model is 

E = C 0 + C j * In (Dose) 



with 



C 0 = 190.45 
Cj = -6.4145 

The submodels can now be assembled and result. in the 
following full model. 



Range = 2368 ‘ ^ g6 + ( 190-6 . u • £n(Dose ))• In (Yield ) (&<2) 

Dose'* 

What the model needs is a function of range and yield 
which yields dose. equation A. 2 cannnot be separated to do 
this, therefore an approximation is offered. The term; 

190 - 6.4 * In (Dose) 

is the slope of the line for a given dose. For doses from 
30 to 10,000 rads, the average slope is 150. The proposal 
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Range (Meters) 



is to use 150 as the slope for every line. Fig. A. 3 has the 
data plotted again as well as Eq. A. 2 and the new proposal. 
For doses of 3,000 rads and below the new proposal seems to 
be an improvement. 




Figure A. 10 Neutron Radiation and Regression. 



Table VIII lists the full model slope and the percent 
error if a slope of 150 is used. Notice than the error is 
negative ar. low doses but is getting larger. Thus, it is 
expected to be in perfect agreement at some dose. Thar dose 
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TABLE VIII 
Slope Approximation 



Dose 


190 - 6. 4* In (Dose) 


Change from 


30 


168 


-12% 


ICO 


160 


-0 6% 


300 


153 


-02% 


516 


150 


00% 


1000 


146 


03% 


3000 


138 


08% 


10000 


131 


13% 



is 518 rads. As with gamma radiation, it is important that 
the model he accurate at radiation levels of 650 rads and 
below. Obviously, this is. 

Above 518 rads, the error is positive and continues to 
increase. At 10,000 rads the error is about the same as at 

30 rads. This seems to be a warning about extrapolation. 

It would be unwise to use this formula for ranges and 
radiation levels outside rhose from which it was derived. 

Using the approximation. Eg. A. 2 becomes 



O o (2 Q h 7 

Range = — ■ ■■ ' 16g6 + 150 • in (.Yield) 

Dose 

which can new be separated. 

The rearranged formula is; 



Dose 



'2368.47 \ 

R-150 • £n (.Yield) ) 



5 .896 



(A. 3) 
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C. THERMAL RADIATION 



The model uses a thermal radiation equation which yields 
range as a function of yield and radiant exposure. Fig. 2.2 
is already in that format with the exception that range is 
in miles. Fig. A. 1 1 is the same figure using a metric range 
scale . 




Figure A. 1 1 Thermal Radiation. 



The previously used method of using submodels will again 
be used. The full model is; 

Range = A + B * In (Yield) 

for each exposure. The coefficients A and 3 are obtained 

for each exposure and are listed in table IX . 

For each exposure, the coefficients are plotted as shewn 
in Figs. A. 1 2 and A. 13. The coefficients can now be 
regressed against the exposure. 
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*1 






TABLE 


IX 


- 




Regression 


Coefficients A and B 






Exposure 


A 


B 






3 


6. 8331 


.5136 






5 


6.5 728 


.505 1 






8 


6.3238 


.5239 






1 2 


6. 1 137 


.5180 






25 


5.7786 


.4926 






50 


5.4593 


.4773 


1 




Fcr 



with 



the first coefficient.. A, the model is 
A = C 0 + C T * In (exposure) 



C 0 = 7.3523 
C = -.4 885 
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Figure A. 13 Regression Coefficient B (Data) . 



The second coefficient could now be regressed against 
exposure, but several models tried failed to yield an 
accurate explanation of behavior of the coefficient as a 
function of exposure. This makes sense since 3 is the slope 
of a given line and all lines appear to be parallel. 
Therefore, an average slope value is used. 

3 = .5051 

The submodels can now be assembled and result in the 
following full model. 



Range = 1561* (Yield) ' 5051 ( A * 4 ) 

v .4835 
(Exposure ) 

Normally, a fig ore would be used to compare the data to 
the regression lines at this point. However, the agreement 
is so good that a graph would lack enough detail to reveal 
any differences. Therefore, Tables X and XI will be offered 
to shew hew well the curve fit is. 
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TABLE X 
Original Data 



pcsure 


1 


2 


Yield 

3 


5 


7 


10 


3 


933 


1289 


1609 


2253 


2594 


2896 


c 


708 


10 14 


1271 


1609 


1931 


2253 


8 


563 


805 


982 


1287 


1496 


1931 


12 


450 


644 


805 


1046 


1255 


1464 


25 


322 


451 


563 


724 


837 


998 


50 


24 1 


322 


386 


499 


612 


708 



TABLE XI 

Regression Results 



Yield 



pcsure 


1 


2 


3 


5 


7 


10 


3 


912 


1295 


1590 


2058 


2439 


2920 


c 


7 1 1 


1009 


1239 


1603 


1900 


2275 


8 


565 


302 


985 


1274 


15 10 


1809 


12 


464 


658 


308 


1045 


1239 


1484 


25 


324 


460 


564 


730 


866 


1037 


5 C 


230 


328 


402 


521 


617 


739 



As cne can clearly see, the 
same caution about extrapolation 
is valid only within the range of 



fit is excellent, but the 
must apply. Equation A. 4 
the data. 
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APPENDIX b 
BOD EL LISTING FILE 



This appendix contains a copy of the listing file for a 
test run of the model. There are three targets, five 
batteries and three companies. Each target has 50 tanks and 
50 trcops uniformly distributed within the RT. Table XII 
lists the target data while Table XIII lists the battery 
data and Table XIV lists the company data. Target types are 
as listed cn page 68. 



TABLE XII 
Target Data 



Target 

1 

2 

3 



Type X 

Coo rdinate 

1 9500.0 

4 8000.0 

5 6000.0 



Y RT 

Coordinate 

0.0 1400.0 

0.0 400.0 

0.0 300.0 
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TABLE XIII 
Battery Data 



Eattery 



Size 


X 

Coordinate 


Y 

Coordinate 


SRC 


0.0 


0.0 


MRC 


1000.0 


0. 0 


SRC 


100 0.0 


500. 0 


aRC 


500.0 


0. 0 


SRC 


1500.0 


250.0 



TABLE XIV 
Company Data 



pany 



Coordinate 



Y 

Coordinate 



2000. C 500.0 
2000.0 0.0 
2000.0 -500.0 



Radius 



200.0 

200.0 

200.0 




100 



04 

ou 

03 

KZ 

6- • 



O O 




2 


M 


2 2 




XM 


O 


3 M 




X 


22 


O M 




*Oj 


XM * 


05 W 




32 


MO * 


• M 




OM 


»W<5 


CJ W 4 




M • 


MM WO 


33 




UM 


0403 


2M X 




X 


Q 


Qh 


% 


NX 


xxx 


4X050 


0 


• 


2 * 


22 3 


X 


W * 


xo-ea 


O .XX 


4 


23 


MM4 


tnO 


O 


XM 


* 4 


O * % 


• 


24 


MW M-t 


O MM 


X 


M4 


oocnx 


3Xt>34 


2 


. 


Q20M 


x cn 


X 


-«M 


X02 


ox . 


M 




43 * 


coz Ou. 




*x 


2 40 CJ 


XX *c 


X 


a. 


X 3053 


MO 




O *M 


.2 


X U X 05 


*M 04 4 
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APPENDIX C 
GLOSSABY 



A. WEAPONS EFFECTS 

Alpha particle: A radioactive emmision identical with the 

helium nucleus, having a mass of four and an electric charge 
of twc positive units. 

Attenuation: A decrease in intensity due to absorption or 

scattering, but not due to geometric reduction as a result 
of range. 

Beta particle: A radioactive emmision identic!® in 

composition to an electron with the possible exception that 
beta particles possess a very high speed. 

Induced Badioactivit y: Radioactivity produced in certain 

elements by the capture of free neutrons. In the area 
immediately below Ground Zero the elements of sodium, 
manganese, aluminum and silicon are are easily induced. 
Neutrcn: A particle of mass one and without any electric 

charge. Neutrons are present in all elements except light 

hydrogen and are required to initiate the detonation process 
in nuclear weapons. 

Bad: Cne of many measures of radiation density or 

destructive power. It represents the absorption of 100 ergs 
of ionizing radiation per gram of absorbing material, such 
as body tissue. 

B. TARGET ANALYSIS 

Functionally impaired: Personnel who exhibit some decreased 

ability to perform their assigned task, but are not 
incapacitated. Characterized by vomiting, diarrhea, nausea, 
lethargy, depression and mental disorientation. 
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Immediate permanent incapacitation (IP) Immediate permanent 
incapacitation is the ;most severe of the radiation criteria 
and causes death in the shortest amount of time. Pers'cnnel 
performing tasks which are physically demanding become 
incapacitated within 5 minutes of exposure. They remain in 
this state for 1 to 2 days at whick time death occurs. 
Immediate transient incapacitation (IT) Personnel receining 
this dose become incapacitated within 5 minutes of exposure 
and remain sc for 30 to 45 minutes. Personnel then recover, 
but are functionally impaired until death which occurs in 4 
to 6 days. 

Latent lethality (LL) Personnel exposed to this their lethal 
dose become functionally impaired within 2 hours of 
exposure. More than half of this group will die within 
several weeks. 

Incapacitated: An individual who performs at 50 % or less 

of his pre-irradiation performance level. 

Radiation exposure status: A unit wide rating relating the 

average accumulated dose to the risk warranted exposure of 
pending missions. 

Minimum separation distance (BSD): The distance friendly 

troops must be from the burst to insure less than 2 % chance 
of damage. 

Light damage: Light damage does not prevent the immediate 

use of an item. Seme repair by the user will usually be 
needed to make the item operational. 

Moderate damage: Moderate damage prevents use of an item 

until extensive repairs are made. This degree of damage is 
normally sufficient to deny use of the equipment. In most 
situations, achievement of this degree of damage will be 
sufficient to support tactical operations. 

Severe damage: Severe damage prevents use of an item 

permanently. There may be some situations such as the 
attack on a bridge where severe damage is the only adequate 
degree of destruction. 
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